## Loading required package: foreach
## 
## Attaching package: 'foreach'
## The following objects are masked from 'package:purrr':
## 
##     accumulate, when
## Loading required package: iterators
## 
## Attaching package: 'data.table'
## The following objects are masked from 'package:lubridate':
## 
##     hour, isoweek, mday, minute, month, quarter, second, wday, week,
##     yday, year
## The following objects are masked from 'package:dplyr':
## 
##     between, first, last
## The following object is masked from 'package:purrr':
## 
##     transpose
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.088587 -0.006528 -0.001660  0.005455  0.060812 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0102102  0.0002809   36.35   <2e-16 ***
## test_set$FCH4 0.1937750  0.0096910   20.00   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01042 on 1800 degrees of freedom
## Multiple R-squared:  0.1818, Adjusted R-squared:  0.1813 
## F-statistic: 399.8 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.214399 -0.004156 -0.000870  0.003302  0.066976 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063025  0.0001128    55.9   <2e-16 ***
## result$FCH4 0.5080646  0.0041981   121.0   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008289 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.6703, Adjusted R-squared:  0.6703 
## F-statistic: 1.465e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.056812 -0.006415 -0.001517  0.005257  0.045742 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0098642  0.0002803   35.19   <2e-16 ***
## test_set$FCH4 0.2669940  0.0108776   24.55   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01016 on 1800 degrees of freedom
## Multiple R-squared:  0.2508, Adjusted R-squared:  0.2504 
## F-statistic: 602.5 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.140559 -0.003996 -0.000804  0.003153  0.091726 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0060996  0.0001033   59.07   <2e-16 ***
## result$FCH4 0.5465918  0.0038445  142.18   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007591 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7373, Adjusted R-squared:  0.7372 
## F-statistic: 2.021e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.046781 -0.006555 -0.001258  0.005046  0.062954 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0100890  0.0002652   38.05   <2e-16 ***
## test_set$FCH4 0.2342332  0.0100435   23.32   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009778 on 1800 degrees of freedom
## Multiple R-squared:  0.2321, Adjusted R-squared:  0.2316 
## F-statistic: 543.9 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.084343 -0.004111 -0.000991  0.003176  0.124377 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0062639  0.0001049   59.72   <2e-16 ***
## result$FCH4 0.5355601  0.0039053  137.14   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007711 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.723,  Adjusted R-squared:  0.723 
## F-statistic: 1.881e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.101141 -0.006816 -0.001937  0.005463  0.065007 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0105040  0.0002833   37.08   <2e-16 ***
## test_set$FCH4 0.2200137  0.0101462   21.68   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01049 on 1800 degrees of freedom
## Multiple R-squared:  0.2071, Adjusted R-squared:  0.2067 
## F-statistic: 470.2 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.222270 -0.004195 -0.000961  0.003341  0.114902 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063836  0.0001097   58.21   <2e-16 ***
## result$FCH4 0.5234734  0.0040831  128.21   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008062 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.6953, Adjusted R-squared:  0.6952 
## F-statistic: 1.644e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.071778 -0.006547 -0.001534  0.005158  0.072791 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0091863  0.0002951   31.13   <2e-16 ***
## test_set$FCH4 0.3050215  0.0112129   27.20   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01072 on 1800 degrees of freedom
## Multiple R-squared:  0.2913, Adjusted R-squared:  0.2909 
## F-statistic:   740 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.136430 -0.003859 -0.000902  0.003060  0.069183 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.005580   0.000107   52.16   <2e-16 ***
## result$FCH4 0.580003   0.003983  145.61   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007865 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7464, Adjusted R-squared:  0.7464 
## F-statistic: 2.12e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.104647 -0.006499 -0.001978  0.004736  0.054491 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0104210  0.0002751   37.88   <2e-16 ***
## test_set$FCH4 0.2350185  0.0104403   22.51   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.0103 on 1800 degrees of freedom
## Multiple R-squared:  0.2197, Adjusted R-squared:  0.2192 
## F-statistic: 506.7 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.223639 -0.004109 -0.001156  0.003031  0.103082 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0064960  0.0001072   60.61   <2e-16 ***
## result$FCH4 0.5328294  0.0039903  133.53   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007879 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7122, Adjusted R-squared:  0.7122 
## F-statistic: 1.783e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.101082 -0.006468 -0.001583  0.004925  0.052405 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0104478  0.0002661   39.26   <2e-16 ***
## test_set$FCH4 0.2210462  0.0099385   22.24   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009968 on 1800 degrees of freedom
## Multiple R-squared:  0.2156, Adjusted R-squared:  0.2151 
## F-statistic: 494.7 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.224632 -0.004161 -0.001050  0.003044  0.115431 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.006478   0.000107   60.52   <2e-16 ***
## result$FCH4 0.530006   0.003985  132.99   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007869 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7106, Adjusted R-squared:  0.7105 
## F-statistic: 1.769e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.065688 -0.006424 -0.001524  0.004955  0.054633 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0101566  0.0002679   37.91   <2e-16 ***
## test_set$FCH4 0.2342388  0.0102115   22.94   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009869 on 1800 degrees of freedom
## Multiple R-squared:  0.2262, Adjusted R-squared:  0.2258 
## F-statistic: 526.2 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.165873 -0.004055 -0.001017  0.003211  0.104060 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0062579  0.0001052   59.51   <2e-16 ***
## result$FCH4 0.5378033  0.0039151  137.37   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007731 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7237, Adjusted R-squared:  0.7237 
## F-statistic: 1.887e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.052091 -0.006942 -0.001918  0.004887  0.075204 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.009605   0.000320   30.02   <2e-16 ***
## test_set$FCH4 0.305081   0.012653   24.11   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.0115 on 1800 degrees of freedom
## Multiple R-squared:  0.2441, Adjusted R-squared:  0.2437 
## F-statistic: 581.4 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.090522 -0.003988 -0.001092  0.002956  0.097211 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0055438  0.0001109   49.98   <2e-16 ***
## result$FCH4 0.5943966  0.0041298  143.93   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008155 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.742,  Adjusted R-squared:  0.7419 
## F-statistic: 2.072e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.094260 -0.006977 -0.001810  0.005738  0.063195 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0107564  0.0002876   37.40   <2e-16 ***
## test_set$FCH4 0.1842998  0.0096624   19.07   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01089 on 1800 degrees of freedom
## Multiple R-squared:  0.1681, Adjusted R-squared:  0.1677 
## F-statistic: 363.8 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.231763 -0.004264 -0.001048  0.003292  0.088957 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063058  0.0001207   52.24   <2e-16 ***
## result$FCH4 0.5282328  0.0044939  117.55   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008873 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.6573, Adjusted R-squared:  0.6572 
## F-statistic: 1.382e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.099896 -0.006180 -0.001648  0.005092  0.063468 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0105304  0.0002713   38.81   <2e-16 ***
## test_set$FCH4 0.2076361  0.0098542   21.07   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01009 on 1800 degrees of freedom
## Multiple R-squared:  0.1979, Adjusted R-squared:  0.1974 
## F-statistic:   444 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.224931 -0.004125 -0.000962  0.003276  0.085707 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0064549  0.0001093   59.08   <2e-16 ***
## result$FCH4 0.5204512  0.0040680  127.94   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008033 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.6944, Adjusted R-squared:  0.6943 
## F-statistic: 1.637e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.066378 -0.006345 -0.001483  0.004949  0.061575 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0102235  0.0002733   37.41   <2e-16 ***
## test_set$FCH4 0.2338459  0.0101743   22.98   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009994 on 1800 degrees of freedom
## Multiple R-squared:  0.2269, Adjusted R-squared:  0.2265 
## F-statistic: 528.3 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.164817 -0.003991 -0.000915  0.003311  0.122674 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0062297  0.0001059   58.81   <2e-16 ***
## result$FCH4 0.5325967  0.0039436  135.05   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007787 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7169, Adjusted R-squared:  0.7168 
## F-statistic: 1.824e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.081565 -0.006991 -0.001949  0.005499  0.096607 
## 
## Coefficients:
##               Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.010179   0.000324   31.41   <2e-16 ***
## test_set$FCH4 0.278310   0.012328   22.58   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01172 on 1800 degrees of freedom
## Multiple R-squared:  0.2207, Adjusted R-squared:  0.2202 
## F-statistic: 509.7 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.174728 -0.004079 -0.001184  0.003038  0.138382 
## 
## Coefficients:
##             Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.005721   0.000114   50.16   <2e-16 ***
## result$FCH4 0.580983   0.004246  136.83   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008384 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7221, Adjusted R-squared:  0.7221 
## F-statistic: 1.872e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.045889 -0.006677 -0.001710  0.004986  0.066562 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0102473  0.0002784   36.80   <2e-16 ***
## test_set$FCH4 0.2612493  0.0109556   23.85   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01022 on 1800 degrees of freedom
## Multiple R-squared:  0.2401, Adjusted R-squared:  0.2396 
## F-statistic: 568.6 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.089230 -0.004130 -0.000944  0.003110  0.124629 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063178  0.0001041   60.69   <2e-16 ***
## result$FCH4 0.5444722  0.0038757  140.48   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007653 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7326, Adjusted R-squared:  0.7325 
## F-statistic: 1.974e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.084042 -0.006504 -0.001639  0.005549  0.038488 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0103947  0.0002558   40.63   <2e-16 ***
## test_set$FCH4 0.1863252  0.0082707   22.53   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.009624 on 1800 degrees of freedom
## Multiple R-squared:  0.2199, Adjusted R-squared:  0.2195 
## F-statistic: 507.5 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.202545 -0.004336 -0.000982  0.003647  0.075656 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0066319  0.0001111   59.68   <2e-16 ***
## result$FCH4 0.4825550  0.0041375  116.63   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.00817 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.6538, Adjusted R-squared:  0.6537 
## F-statistic: 1.36e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.054633 -0.006581 -0.001505  0.004878  0.062886 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0101317  0.0002792   36.28   <2e-16 ***
## test_set$FCH4 0.2488734  0.0103963   23.94   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01011 on 1800 degrees of freedom
## Multiple R-squared:  0.2415, Adjusted R-squared:  0.2411 
## F-statistic: 573.1 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.142137 -0.003974 -0.000927  0.003076  0.110841 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0060593  0.0001046   57.94   <2e-16 ***
## result$FCH4 0.5414592  0.0038935  139.07   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007688 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7286, Adjusted R-squared:  0.7286 
## F-statistic: 1.934e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.060266 -0.006622 -0.001702  0.005191  0.053183 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0105263  0.0002763   38.09   <2e-16 ***
## test_set$FCH4 0.2795353  0.0104466   26.76   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.0101 on 1800 degrees of freedom
## Multiple R-squared:  0.2846, Adjusted R-squared:  0.2842 
## F-statistic:   716 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.138062 -0.004134 -0.001171  0.003135  0.087972 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063390  0.0001023   61.95   <2e-16 ***
## result$FCH4 0.5414730  0.0038100  142.12   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007523 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7371, Adjusted R-squared:  0.7371 
## F-statistic: 2.02e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.108971 -0.006238 -0.001503  0.004592  0.062790 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0104501  0.0002724   38.37   <2e-16 ***
## test_set$FCH4 0.2390721  0.0101679   23.51   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01005 on 1800 degrees of freedom
## Multiple R-squared:  0.235,  Adjusted R-squared:  0.2345 
## F-statistic: 552.8 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.226219 -0.004132 -0.001025  0.003123  0.123061 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063691  0.0001059   60.16   <2e-16 ***
## result$FCH4 0.5330332  0.0039418  135.22   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007783 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7174, Adjusted R-squared:  0.7173 
## F-statistic: 1.829e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.094774 -0.006536 -0.001607  0.004900  0.041756 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0101341  0.0002719   37.27   <2e-16 ***
## test_set$FCH4 0.2004374  0.0093309   21.48   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01016 on 1800 degrees of freedom
## Multiple R-squared:  0.204,  Adjusted R-squared:  0.2036 
## F-statistic: 461.4 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.218365 -0.004080 -0.001016  0.003417  0.082211 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0063685  0.0001111   57.31   <2e-16 ***
## result$FCH4 0.5090022  0.0041372  123.03   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.008169 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.6775, Adjusted R-squared:  0.6775 
## F-statistic: 1.514e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Joining with `by = join_by(FCH4, Ustar, NEE, LE, H, Rg, Tair, RH, VPD, DoY,
## Ts_1_5cm, Ts_1_10cm, Ts_1_20cm, Ts_1_50cm, Cond, WL, TW, Ts_2_5cm, Ts_2_10cm,
## Ts_2_20cm, Ts_2_50cm, s, c)`

## [1] "RF Summary of test data:"
## 
## Call:
## lm(formula = test_set$FCH4_rf ~ test_set$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.068610 -0.006543 -0.001962  0.004819  0.084404 
## 
## Coefficients:
##                Estimate Std. Error t value Pr(>|t|)    
## (Intercept)   0.0089865  0.0003133   28.68   <2e-16 ***
## test_set$FCH4 0.3193849  0.0122425   26.09   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.01129 on 1800 degrees of freedom
## Multiple R-squared:  0.2744, Adjusted R-squared:  0.274 
## F-statistic: 680.6 on 1 and 1800 DF,  p-value: < 2.2e-16

## [1] "RF Summary of whole data:"
## 
## Call:
## lm(formula = result$FCH4_RF_model ~ result$FCH4)
## 
## Residuals:
##       Min        1Q    Median        3Q       Max 
## -0.096717 -0.003831 -0.000957  0.003028  0.098357 
## 
## Coefficients:
##              Estimate Std. Error t value Pr(>|t|)    
## (Intercept) 0.0054318  0.0001081   50.27   <2e-16 ***
## result$FCH4 0.5919832  0.0040230  147.15   <2e-16 ***
## ---
## Signif. codes:  0 '***' 0.001 '**' 0.01 '*' 0.05 '.' 0.1 ' ' 1
## 
## Residual standard error: 0.007944 on 7204 degrees of freedom
##   (8631 observations deleted due to missingness)
## Multiple R-squared:  0.7504, Adjusted R-squared:  0.7503 
## F-statistic: 2.165e+04 on 1 and 7204 DF,  p-value: < 2.2e-16
## Warning: `as.tibble()` was deprecated in tibble 2.0.0.
## ℹ Please use `as_tibble()` instead.
## ℹ The signature and semantics have changed, see `?as_tibble`.
## This warning is displayed once every 8 hours.
## Call `lifecycle::last_lifecycle_warnings()` to see where this warning was
## generated.
##             StartDate             EndDate AvailNum   EffNum ObsNum AnnNEE_gCO2
## 1 2022-06-01 15:00:00 2023-04-27 13:00:00    15837 1477.062   9930    29.82144
##   UncNEE_Rand_gCO2 UncNEE_uStar_gCO2 UncNEE_Comb_gCO2 AnnNEE_gC UncNEE_Rand_gC
## 1         63.53007          32.61282         71.41194  8.138726       17.33832
##   UncNEE_uStar_gC UncNEE_Comb_gC
## 1        8.900535       19.48941
##             StartDate             EndDate AvailNum   EffNum ObsNum AnnFCH4_gCH4
## 1 2022-06-01 15:00:00 2023-04-27 13:00:00    15837 2238.738   7034        6.101
##   UncFCH4_Rand_gCH4 UncFCH4_uStar_gCH4 UncFCH4_Comb_gCH4 AnnFCH4_gC
## 1         0.2003837          0.1353286         0.2418004   4.568523
##   UncFCH4_Rand_gC UncFCH4_uStar_gC UncFCH4_Comb_gC
## 1       0.1500504        0.1013362       0.1810639